Effect of atomic electrons on 7.6 eV nuclear transition in 229 Th 3+ 



S. G. Porsev 1,2 and V. V. Flambaum 1 



o 
o 

> 

o 



43 . 

i 



03 

1*3 1 

o ■ 

C/5 ■ 
>Y 

43 : 

Ph. 



> 

o 

On 

m 



OS 

o 




School of Physics, University of New South Wales, Sydney, NSW 2052, Australia and 

2 Petersburg Nuclear Physics Institute, Gatchina, Leningrad district, 188300, Russia 

(Dated: November 19, 2009) 

We have considered an effect of atomic electrons due to the electronic bridge process on the 
nuclear 229m Th - 229 9Th transition in 229 Th 3+ . Based on a recent experimental result we assumed 
the energy difference between the isomeric and the ground nuclear states to be equal to 7.6 eV. We 
have calculated the ratios of the electronic bridge process probability (Feb) to the probability of the 
nuclear radiative transition (Tn) for the electronic 5/5/2 — * 6^3/2, 6^5/2, 7s and the 7s — > 7pi/2, 7p3/2 
transitions and found Y-^b./'^n ~ 0.01 -=-0.1 for the former and Feb/IV ~ 20 for the latter. 

PACS numbers: 31.15.Ar, 23.20.Lv, 27.90+b 



I. INTRODUCTION 

The 229 Th nucleus is unique in a sense that the en- 
ergy splitting of the ground state doublet is only several 
eV. Though a prediction of existence of so low-lying level 
has been made more than thirty years ago [l|, the defi- 
nite value of energy of the isomeric state 229rn Th is not 
known so far. In 1990 Reich and Helmer 0] measured 
this excitation energy (ujjf) to be 3.5 ± 1.0 eV. In Ref. [3[ 
it was obtained_5.5± 1.0 eV. Finally, a recent experiment 
of Beck et al. [H has given even larger value (with least 
error) ujn = 7.6 ± 0.5 eV. 

As to the lifetime of the 229m Th, measurements per- 
formed by different experimental groups led to different 
values. The results differ from each other by several or- 
ders of magnitude, changing from a few minutes [5J to 
many hours Q. Hence, new experimental and theoreti- 
cal investigations are required. 

A special interest to the nuclear transition from the iso- 
meric state to the ground state is motivated by a possibil- 
ity to build a superprecise nuclear clock 0] and very high 
sensitivity to the effects of possible temporal variation of 
the fundamental constants including the fine structure 
constant a, strong interaction and quark mass Q. 

Laser cooling of the 232 Th 3+ ion was recently reported 
by Campbell et al. in their paper Q. This was the first 
time when a multiply charged ion has been laser cooled. 
As a next step this experimental group plans to inves- 
tigate the nuclear transition between the isomeric and 
the ground state in a trapped, cold 229 Th 3+ ion. Moti- 
vated by this experimental progress we have considered 
229 Th 3+ ion and calculated the transition probability of 
the 229 Th nucleus from its lowest-energy isomeric state 
229m Th to the ground state 229 9Th due to the electronic 
bridge (EB) process. 

Our calculations, based on the value of ojn = 7.6 eV, 
showed that if the electrons are in their ground state 
the ratio of the probability of the EB process, Teb, 
to the probability of the nuclear radiative Ml tran- 
sition, rjv(Ml), is of the order of (a few)xl0^ 2 . If 
the valence electron is in the metastable 7s state then 



FIG. 1: The single and double solid lines relate to the elec- 
tronic and the nuclear transitions, correspondingly. The 
dashed line is the photon line. 



r EB /r w (Mi) ~ 20. 

The paper is organized as follows. In Section |n] we 
present the general formalism describing the EB process. 
Section Mil is devoted to the method of calculation of 
the properties of Th 3+ . In Section IIVI we discuss the 
results of calculations and Section fVl contains concluding 
remarks. If not stated otherwise the atomic units (H = 
|e| = m e = 1 and the speed of light c = 137) are used. 



II. GENERAL FORMALISM 



Configurations mixing between combined 
electron-nucleus states 



The 7.6 eV transition in 229 Th is the Ml transition 
with the amplitude of a fraction of the nuclear magneton 
fiN. An amplitude of an allowed electric dipole transition 
of the valence electron, ~ 1 au, is 10 6 times larger. If 
there is an electron excited state close to the energy of the 
nuclear excitation, an energy transfer from the nuclear 
excited state to the electron excited state accompanied 
by the electron electric dipole transition to a lower state, 
may radically decrease the lifetime of the nuclear isomeric 
state. Even if there is no an electron state very close 
to the nuclear excited state, the electron bridge process 
produces significant effect. 

The EB process can be represented by two Feynman 
diagrams in Fig. [TJ In the following we assume that the 
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initial i and the final / electronic states are of opposite 
parity and a real photon which is emitted or absorbed 
is the electric dipole photon. The probability of the EB 
process in this case is much larger than in the case when 
the i and the / states are of the same parity. 

Therefore, the EB process can be effectively treated as 
the electric dipole i — > f transition of the electron accom- 
panied by the nuclear transition from its isomeric state 
to the ground state. Denoting by Deb the amplitude of 
this "generalized" electric dipole transition and assuming 
that the initial and the final states are fixed, we obtain 
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where the indices i, (n, fc), and / denote initial, interme- 
diate, and final electronic states, correspondingly; and 
the indices g and m denote the ground state and the iso- 
meric state of the nucleus. Si are the atomic energies, 
E m i g \ are the nuclear energies of the isomeric (ground) 
states, and Ti are the widths of the intermediate states 
which may be neglected in the case of 229 Th. The op- 
erator D = — r is the electron electric dipole moment 
operator and H mt is the hyperfine coupling Hamiltonian, 
which may be represented as a sum over multipole nu- 
clear moments M.\ of rank K combined with the even- 
parity electronic coupling operators Tkx of the same rank 
as 
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Neglecting the hyperfine splitting of levels, we can rep- 
resent the total wave function as a product of the nuclear 
wave function and the electronic wave function. For in- 
stance, \g,n) = \g)\n) = \I g M g )\j n J n m n ), where I g is 
the nuclear spin, M g is the projection of the nuclear 
spin; J n is the electron total angular momentum, m n 
is its projection, and j n encapsulates all other electronic 
quantum numbers. Taking into account Eq. ([2]) we can 
rewrite Eq. |T]) as 
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Eg. 



where uj a b = e a — £b and ojn = E„ 

Thus, we need to carry out the atomic calculation 
which is similar to that for a forbidden El transition 
opened by the hyperfine interaction (see, e.g., [HI)- The 
only difference is that the matrix element (ME) of the 
nuclear moment (<7|.M^-|m) here is non-diagonal (there 
is also a few per cent correction due to variation of the 
electron magnetic field inside the nucleus). Note that 



the conventionally defined nuclear moments are related 
to the tensors M K as fi = (IM I = I\M°\IMj = I) and 
Q = 2{IMi = I\M%\IMi = I). 

The probability Teb of the electric dipole transition 
determined by its amplitude Deb is given by a simple 
formula (see, e.g., (i"Tj |) 
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where u> is the real photon frequency determined from 
the low of conservation of energy as lo = e% — sj + lun- 

If we average over the initial projections of the elec- 
tronic and the nuclear total angular momenta m.j and 
M m and summing over the final projections m/ and M g , 
Eq. (UJ) is transformed to 
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Substituting Eq. (O to Eq. ([5]), applying the Wigner- 
Eckart theorem and performing the summation over all 
magnetic quantum numbers of the initial, intermediate 
and final states, we can reduce Eq. ([5]) to the form Teb = 

J2k Tee?' wriere Tees' 1 can ^ e represented by 
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The terms G[ K ^ and G^ characterize the contributions 
of the first and second diagrams in Fig. [1] while the "in- 
terference" of two these diagrams is given by G^'. 

It is worth noting that Eq. ([6]) is valid in a general 
case because deriving it we did not make any approxima- 
tions. In particular, we did not suppose that there is an 
electronic transition whose frequency is close to the the 
nuclear transition frequency lun- In systems where such 
a "resonance" transition exists, the expression for Teb 
can be significantly simplified. 



B. 



Derivation of the coefficients @mi and (3e2 



Since the frequency of the nuclear transition from the 
isomeric state to the ground state of 229 Th is very small, 
in the following we will take into consideration only first 
two terms in Eq. ([6]), involving the nuclear magnetic- 
dipole (K — 1) and electric-quadrupole (K = 2) mo- 
ments. Another consequence of the smallness of the 
nuclear transition frequency is that the probability of 

the m g transition Tn(E2) is strongly suppressed 

in comparison to the probability of m ^—f q transition 
T N (M1). 

The probability T jv (tK , m — > g) of the tK transition 
(where r denotes M or E) in the 229 Th nucleus can be 
written in a form used in the nuclear physics as (see, 
e.g., 0) 



Tn{tK, m — > g) 
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[(2i^ + l)!!] 2 K 
ujn/c and the reduced probability of the nu- 



Here 

clear m — > g transition B{tK,tti 
terms of the operator Mk, reads as 



g), expressed in 



B(tK, m — > g) = 
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Using Eq. (fT0|) we find for this transition 
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The theoretical value of B(Ml,m — > g) 
in [l| 

B(M1, m -> g) w 0.086 (j? n . 



9)' 

was obtained 
(13) 



To the best of our knowledge the accurate value of 
B(E2,m — > g) is unknown. An estimate of this quantity 
is found in Ref. Il4l. where Strizhov and Tkalya, refer- 
ring to the paper [l5| , cite the value of several Weisskopf 
units (W.u.) for B(E2,m — > g). 

The definition of 1 W.u. for the E2 transition from a 
nuclear excited state to the ground state (in usual units) 
is 



E2 : 1 W.u. = 5.940 x 10~ 6 ,4 4/3 (e • barn) 2 



where e is the electron charge and A is the number of 
nucleons in the nucleus. 

Using Eqs. dH]), |T5]), and (fj?]) we arrive at the esti- 
mate 
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(15) 



An accurate calculation of the probability of the nuclear 
E2 transition is beyond the topic of this work. In the 
following we rely on the estimate given by Eq. (|15|) and 
concentrate our efforts on the computation of the ratios 

0mi = r^/rjv(Ml) and /3 E 2 = T ( ^l/T N {E2), 

where T^g 2 ^ are given by Eq. ([6]) and Tjv(MI) and 
T N {E2) can be found from Eq. (fT0| . 
Using these equations we obtain 
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As follows from the estimate Eq. (115)) , the probability 
of the nuclear radiative E2 transition from the isomeric 
state to the ground state in 229 Th is completely negligible 
in comparison with the probability of the Ml transition. 
Based on this estimate one can expect that the electronic 
part of the EB process mainly contributing to Teb can 
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f, while the channel 



(14) 



be represented as i 

n — -> / can be neglected. 

As we will demonstrate below this assumption is valid 
for 229 Th 3+ in spite of that (3e2 is many orders of mag- 
nitude larger than (3mi- The physical meaning of this 
is as follows. It is known that a neutral atom is not 
affected by an external electric field. It means that an 
effective electric field acting on the nucleus is equal to 
zero because the electrons completely screen the exter- 
nal electric field. Respectively, gradient of electrostatic 
potential created by the electrons at the nucleus is very 
large. For a static case (in our consideration it corre- 
sponds to lun = 0) a similar phenomenon was investi- 
gated in [l|| [ItJ where magnetic-dipole shielding factors 
and electric-quadrupole antishielding factors were calcu- 
lated for a number of atoms and ions. For instance, for 
such a heavy atom as Hg, the latter was shown to be four 
orders of magnitude larger than the former. 

Note also that the probability of, so called, "elastic" 
process (when the final state is the same as the initial 
state) is much smaller, since instead of El transitions 
we have to consider Ml (or E2) transitions. But the 



probability of an allowed Ml transition is four orders of 
magnitude smaller than the probability of an allowed El 
transition. 

The triply ionized thorium 229 Th 3+ is an univalent ion. 
Respectively, the total electronic angular momentum as 
well as other quantum numbers coincide with the quan- 
tum numbers of the valence electron. The expressions 
for the single-electron operators T\ and Ti and for the 
MEs of the operators D, T\ and X2 are presented in the 
Appendix [A] 

III. METHOD OF CALCULATION 
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FIG. 2: The single and double solid lines relate to the elec- 
tronic and the nuclear transitions, correspondingly. The 
dashed line is the photon line. The wavy line is the Coulomb 
interaction between electrons. 



At the first stage we have solved Dirac-Hartree-Fock 
(DHF) equations [l8| in V N ~ X approximation. It means 
that the DHF equations were solved self-consistently 
for the core electrons. After that we determined va- 
lence orbitals for several low-lying states from the frozen- 
core DHF equations. The virtual orbitals were deter- 
mined with the help of a recurrent procedure [l9j . One- 
electron basis set of the following size was constructed: 
1 - 20s, 2 - 20p, 3 - 20d, 4 - 25/, 5 - I85. 

To find wave functions needed for calculation of ft/i 
and Pe2 we applied a relativistic many-body method 
initially suggested in Refs. 0, [HJ and subsequently 
developed in [2^, HU. In this method one determines 
wave functions from solution of the effective many-body 
Schrodinger equation 

H eS (E n )\y n )=E n \y n ), (18) 

with the effective Hamiltonian defined as 



H eff (E) = H FC + £(£) 



(19) 



Here H-pc is the frozen-core DHF Hamiltonian and self- 
energy operator £ is the energy-dependent correction, 
involving core excitations, which recovers second order 
of perturbation theory in residual Coulomb interaction 
and additionally accounts for certain classes of many- 
body diagrams in all orders of perturbation theory. We 
will refer to this approach as the DHF+S formalism. 

Together with the effective Hamiltonian H c g we in- 
troduce effective ( "dressed" ) electric-dipole operator D e g 
and operators (Tk) c s acting in the model space of valence 
electrons. These operators were obtained within the 
relativistic random-phase approximation (RPA) [Tt], [HJ 
which describes a shielding of the externally applied elec- 
tric field by the core electrons. The RPA sequence of 
diagrams was summed to all orders of the perturbation 
theory. 

A representative diagram illustrating a contribution of 
the RPA corrections in the first order is shown in Fig. [2] 
As we will show below in certain cases including the RPA 
corrections is very important because it changes Teb by 
orders of magnitude. 

With the wave functions obtained from Eq. (TTS)) . the 



the Sternheimer [24( or Dalgarno-Lewis [25| method im- 
plemented in the DHF+RPA+S framework. 



For instance, the expression for G 2 
can be rewritten as 



(K) 



given by Eq. 
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where an intermediate-state wave function \5^) can be 
found from the inhomogeneous equation 
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off 
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quantities G\ , G^\ and G 2 can be computed with 



and then \5ip, J n ) is obtained by projecting the wave func- 
tion I (5-0) to the state with the definite value of J n . Sim- 
ilarly we can derive the expressions for and G^P . 

Only excitations of the valence electron to higher vir- 
tual orbitals are included in the intermediate-state wave 
function |<5-0) due to the presence of H e s in Eq. (|2ip . 

Additional contributions to G^, G^, and G% come 
from particle-hole excitations of the core. The role of 
these contributions will be discussed more detailed in the 
next section. 

Since Th 3+ is an univalent element, the quantities 
G^ "^j G^2 \ and G 2 can be obtained by another 
method. We can directly sum over all intermediate 
states using the single-electron wave functions found at 
the stage of constructing the basis set. An accuracy of 
this approach is comparable to the accuracy of the more 
refined method of solving the inhomogeneous equation. 
The reason is that, despite a non-resonant character of 
the EB process in 229 Th 3+ for lon — 7.6 eV there are only 
a few intermediate states in Eqs. (|7|), (JSJ), and (O (whose 
denominators are small) that give a dominant contribu- 
tion to Teb- 

We would like to stress that in the sums over the in- 
termediate states in Eqs. (O, ([8]), and ([9]) the states 
{liJi) = (jkJn) or (7/J/) = (7fc J n ) are allowed. This 
is due to that the ME of the nuclear moment (g\Ai\\m) 
is non-diagonal and, correspondingly, the diagonal MEs 
of the electronic operator Tk should be included into con- 
sideration. Note that the diagonal MEs of the operators 
Tk are large and the inclusion of these contributions to 
Teb significantly affects the final value of the latter. 



5 



TABLE I: The low-lying energy levels (in cm -1 ) in the DHF 
and the DHF+E approximations are presented. The theoret- 
ical values are compared with the experimental data. 





DHF 


DHF+E 


PrVTlPT'lTVIPTli - a 
UAUCi 1111C11L 


OO3/2 




0/5/2 






DU5/2 










0/5/2 


oiyu 


003/2 


yoyi 


y±yo 


OJ7/2 


sal 7 


605/2 


14ooo 


1/1/1 Qd 


'•=1/2 


1 1 ^1 Q 


7Sl/2 


21321 


9*31 *31 

ZrO-LO-L 


7P1/2 


46702 


7P1/2 


59436 


60239 


7P3/2 


58225 


7P3/2 


72677 


73056 


8si /2 


102595 


8fll/2 


120085 


119622 


7d 3 /2 


103148 


7^3/2 


120898 


119685 


7d 5 /2 


104763 


7^5/2 


122657 


121427 


6/5/2 


111874 


6/5/2 


128734 


127262 


6/7/2 


112316 


6/7/2 


129227 


127815 


8P1/2 


117185 


8P1/2 


135144 


134517 


8P3/2 


122194 


8P3/2 


140558 


139871 


9S1/2 


142328 


9Sl/ 2 


161481 


160728 






"Reference 


271; 





6 The removal energy of the 5/5/2 state was found to be 
equal to 0.9414 au on the DHF stage and 1.0584 au on the 
(DHF+E) stage. The experimental value is 1.0588 au. 



IV. RESULTS AND DISCUSSION 

To check the quality of the constructed wave functions 
we have calculated the energy levels for a number of low- 
lying states and compared them with the experimental 
data. Some details regarding the energy levels computa- 
tion can be found in our recent paper (2d| . We present 
in Table [J the results obtained on the stage of pure DHF 
approximation and in the frame of DHF+E formalism. 

As seen from Table fl] on the stage of the DHF ap- 
proximation the order of the low- lying levels is incorrect. 
For instance, the 6(i 3 /2 state lays deeper than the 5/5/2 
state. An agreement between theoretical and experimen- 
tal energy levels is rather poor. The inclusion of the 
core- valence correlations restores the correct order of the 
states and significantly improves the agreement with the 
experimental energy levels. Nevertheless in certain cases 
(e.g., for the 7s state) the energy levels were reproduced 
not very accurately. For this reason in the following cal- 
culation of Pmi and [3e2 we used the experimental ener- 
gies for the low-lying states. 

We have carried out calculations of the coefficients /3mi 
and @E2 for ujn — 7.6 eV considering the ground state 
5/5/2 and the metastable state 7s as the initial state i. As 
follows from the discussion above the final states should 
be of opposite parity in comparison to the initial states. 
Respectively, we considered 6d 3 / 2 , 6d 5 / 2 and 7s states to 
be the final states when the initial state was 5/5/2. The 
7pi/2 and the 7p 3 / 2 states were the final states when the 
initial state was 7s. 

In Table HT1 we present the results obtained 1) on the 
stage of pure DHF approximation, 2) in the DHF + RPA 
approximation, and 3) in the frame of DHF + RPA + 



TABLE II: The coefficients /3a/i obtained for certain i — ► / 
transitions for u N = 7.6 eV in the DHF, the DHF+RPA, 
and the DHF+RPA+E (denoted as +E) approximations are 
presented. The coefficients (3e2 are given in the DHF+RPA 
approximation. 



i 


/ 


DHF 


(3 Ml 

DHF+RPA 


+E 


(3E2 

DHF+RPA 


5/5/2 


603/2 


0.015 


0.0037 


0.015 


2 xlO 8 


6d 5 / 2 


0.0015 


0.051 


0.060 


5 xlO 7 




7Si/ 2 


2xl0~ 9 


0.032 


0.037 


5 xlO 7 


7si/ 2 


Tpi/2 


18 


19 


23 


7 xlO 8 




7f>3/2 


4.1 


4.4 


5.3 


1 xlO 8 



S formalism. Including the RPA corrections is formally 
reduced to replacement of the "bare" operators by the 
"dressed" operators. In particular, solving the inhomo- 
geneous equation we have to replace the operators Tr- in 
Eq. HOI and D in Eq. (HJ) by (T K ) eS and D cS , corre- 
spondingly. 

As seen from the table in certain cases the inclusion 
of the RPA corrections increases the probability of the 
EB process by several orders of magnitude. It happens, 
for example, for the 5/5/2 — ► 7s transition. The channel 

5/5/2 — ^ n 7s turns out strongly enhanced because 
the "dressed" MEs (5/5/ 2 ||(7i) e ff||^} are much larger in 
absolute value than the "bare" MEs (5/5/ 2 ||7i||n). In- 
deed, we have to consider the intermediate states n that 
admit the El transitions n — > 7s. But for such n the 
"bare" MEs |(5/ g / 2 ||7i||n)| are very small. 

The coefficient (3 mi is rather small for the 5/5/2 — > 
6d 3 / 2 transition and the RPA and the E corrections 

change its value significantly. The reason is that 
G$, and G { 2 1] are comparable in their magnitudes but 

G^' is negative. In the DHF+RPA approximation it 
leads to a large cancellation between these terms. 

When we consider the 7s state as the initial state, the 
main channel of the process is 7s — 1 -+ 8s 7px/2,3/2- 
Respectively, the first diagram in Fig. [T] (the term G^) 

gives the main contribution to Teb while and G 2 ^ 
only slightly correct this value. As is seen (3mi(7s — > 7pj) 
are 2-3 orders of magnitude larger than /3mi(5/ 5 / 2 — > 
6dj;7s). This is due to the large value of the ME 
(7s||7i||8s). 

As is seen from Table [H] the inclusion of the core- 
valence correlations changes the values of /3mi at the 
level of 20% for all considered transitions except the 
5/5/2 ~~ * 6g? 3 /2 transition. These corrections are not too 
large because the core orbitals lay rather deep. In partic- 
ular the single-electron energy of the external core 6p 3 / 2 
orbital is —2.1 au. For the same reason the contribution 
to Teb from the core electrons excitations is small. It is 
at the level of few per cent. 

We also present in Table |TT] the coefficients (3e2 ob- 
tained for the 5/5/2 — > 6dj;7s and the 7s — > 7pj tran- 
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TABLE III: The probabilities F^g (in sec" 1 ) obtained for cer- 
tain i — > / transitions for cjjv = 7.6 eV in the DHF+RPA+E 
approximation are presented. 



5hJ~2 9.9 xlO-° 

64/2 4.0 xlO" 5 

7s 1/2 2.4 xlO" 5 

7s 1/2 7p V a 1.5 xlCT 2 

7p 3 /2 3.5 xlO" 3 



contribution of the £ — ^ n — > f channel to Feb is negli- 
gible. It is worth noting that this statement is correct for 
all considered transitions in spite of the resonant charac- 
ter of the 5/5/2 — ^> 7pi/2 — — > 6d 3 / 2 transition. 

We would like to thank J. Berengut and A. Kuzmich 
for stimulating discussion. This work was supported by 
Australian Research Council. The work of S.G.P. was 
supported in part by the Russian Foundation for Basic 
Research under Grants No. 07-02-00210-a and No. 08- 
02-00460-a. 



sitions in the DHF+RPA approximation. We restricted 
ourselves by this simple approximation because these val- 
ues are given mostly for reference and an order of mag- 
nitude estimate of these quantities is sufficient. 

As we have already mentioned above the coefficients 
(3e2 are many orders of magnitude larger than the coeffi- 
cients Pmi found for the same transitions. In particular, 
for the 5/5/2 — > 6rf 3 / 2 transition 0E2/0M1 ~ 10 10 . It 
is not surprising if we take into account the small value 
of in the denominator of Eq. (fT7|) and the resonant 

character of the 5/ 5 / 2 — ^> 7pi/2 6d 3 / 2 transition 
because the frequency of the 5/5/2 - 7pi/2 transition 
UJ 7pi/2, 5/5/2 ~ 7.5 eV is very close to u>n = 7.6 eV. 

In spite of that the main contribution to Feb comes 

from the i n -^i> / channel. As it follows from the 
results listed in Table [TT1 and Eq. (|15[) . we can neglect 

the contribution to Teb coming from the i n / 



(i) 
eb- 



channel and put Teb ~ T 

Using Eqs. CGI]) and {BJ we find T N (M1) « 6.6 x 
10 -4 sec" 1 at cjjv = 7.6 eV and, correspondingly, 



r 



EB 



,(1) 
EB 



6.6 x 10 Pmi sec 



The numerical results obtained for T^g with use of the 
equation written above are listed in Table IIIII 



(22) 



CONCLUSION 



In conclusion, we have calculated the ratios of the prob- 
abilities r E g and Tgg to the probabilities of the nuclear 
radiative Ml and E2 transitions, Pmx an d Pe2- We 
found that if the valence electron is in the ground state 
the coefficients (3mi are rather small for all considered 
transitions. If the valence electron is in the metastable 
7s state the coefficients (3mi are 2-3 orders of magnitude 
larger and r EB /IV(Ml) - 20. 

The spectrum of Th 3+ is not too dense. As a result 

for the i — i+ n / transitions considered in this work 
there are no electronic transitions which would be at res- 
onance with the nuclear transition at ojm = 7.6 eV. 

We have found the coefficients (3e2 to be many orders 
of magnitude larger than (3mi, but based upon the es- 
timate r Ar (Ml)/r7v(-E , 2) ~ 10 11 one can state that the 



APPENDIX A 

The expressions for the single-electron operators T± 
and Ti can be written as 

, , -isfla. ■ df) (n) , . 

Tix{v) = — (Al) 



and 



_ -C 2 A(n) 



(A2) 



where n = r/r and is a normalized vector spherical 
harmonic defined by (see, e.g., [28| l 



y/KW+T) 



C K \(n). 



(A3) 



Here L is the orbital angular momentum operator and 
Ckx is a spherical harmonic given by 



4tt 



2K + 1 



Y K \(n). 



(A4) 



To calculate the MEs of the operators D, T\, and T2 
we define the one-electron wave function \a) = V'a(r) as 
follows 

^ (r) = if % r )& m - {n) , (A5) 

where n a = (l a - j a )(2j a + 1). 
Using the ME («6||Cjr||K„) : 



(h\\Ck\K) = (-l)^ +1 /V(2ja + l)(2j b + l) 

jb ja K 



x 1 1/2 1/2 I t(!b + la + K), 



where 



1, if x is even 
0, if x is odd, 



we can write the reduced ME for the electric dipole op- 
erator D in the following form 

(n b K b \\D\\n a K a ) = (n b K b \ \ - r||n a K a ) = 
-(Kb I 



i||Ci||K a ) / {P b Q a +PaQb}rdr, (A6) 
Jo 



7 



where rij is the principal quantum number. 

The reduced ME for the magnetic dipole operator T\ 
is represented by 

(n b K b \\Ti\\n a K a ) = (-Kf,||Ci||K a ) 

f°° 1 
x(n b + K a ) / {P b Qa+PaQb} — dr. (A7) 
Jo r " 

Rewriting the angular part of Eq. (|A7|) in a more simple 
form we arrive at 

(n b K b \\T 1 \\n a K a ) =ah + l a )(-l) ja+la - 1/2 

c ■ r°° 1 

x / { PbQa +p a Q b } — dr, (A8) 

2 Jo r 



where c jajh . ( ^ ± ,+ D/O'^+l). * / * 

and j min = min (j a ,j b ). 

The reduced ME for the electric quadrupole operator 
X2 is given by 



(n b K b \\T 2 \\n a K a ) = -(n b \\C 2 \\Ka} 

x / {P a Pb +Q a Q b }^ dr. (A9) 
Jo r 
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